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Aims. An evolution of the electron distribution function in the beam-plasma system with the return current is com- 
puted numerically for different parameters. The X-ray bremsstrahlung corresponding to such an electron distribution 
is calculated and the directivity of the X-ray emission is studied. 

Methods. For computations of the electron distribution functions we used a 3-D particle-in-cell electromagnetic code. 
The directivity of the X-ray emission was calculated using the angle-dependent electron-ion bremsstrahlung cross- 
section. 

Results. It was found that the resulting electron distribution function depends on the magnetic field assumed along 
the electron beam propagation direction. For small magnetic fields the electron distribution function becomes broad 
in the direction perpendicular to the beam propagation due to the Weibel instability and the return current is formed 
by the electrons in a broad and shifted bulk of the distribution. On the other hand, for stronger magnetic fields the 
distribution is more extended in the beam-propagation direction and the return current is formed by the electrons in 
the extended distribution tail. In all cases, the anisotropy of the electron distribution decreases rapidly due to fast 
collisionless processes. However, the magnetic field reduces this anisotropy decrease. The X-ray directivity shows the 
same trend and it is always closer to the isotropic case than that in a simple beaming model. 
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1. Introduction 



It is commonly believed that the hard X-ray emission in 
solar flares is produced by the bremsstrahlung process of 
energetic electrons in dense layers of the solar atmosphere 
(Brown 1971; Tandberg-Hanssen & Emslie 1988). 

It is also known that up to now this scenario has sev- 
eral unresolved drawbacks as summarized in the paper 
by Brown et al. (1990). For example, the bremsstrahlung 
mechanism generating the hard X-ray bursts is of a very 
low efficiency and therefore huge electron beam fluxes Ep 
= 10^ - 10^^ ergs s~^ cm~^ are required for an explanation 
of the observed X-ray fluxes (Hoyng et al. 1978). It means 
that at the acceleration site in the low corona with a rel- 
atively low density (rie ~ 10^ cm~^), a substantial part of 
all plasma electrons needs to be accelerated. Furthermore, 
these electron beams represent huge electric currents that 
have to be neutralized by the return currents. The return 
current is a natural part of any beam-plasma system (van 
den Oord 1990). 

The beam-plasma interaction has been studied for a 
long time, starting with the paper by Bohm & Gross (1949). 
While the flrst 1-D models considered the electrostatic as- 
pects of this interaction (two-stream instability, generation 
of Langmuir waves, and quasi-linear relaxation of the beam, 
see e.g. Melrose 1980, Birdsall & Langdon 1985, Benz 1993, 
Karlicky 1997 and the references therein), new 3-D stud- 
ies include the return current and electromagnetic effects 
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which lead to many further instabilities (Weibel, filamenta- 
tion, oblique, Bell, Buneman, and so on, see Karlicky 2009, 
Bret 2009). (Remark: The Weibel instability in the sense 
used here and in the paper by Nishikawa et al. (2008) is 
also known as the filamentation instability (Bret 2009).) 
To cover all these processes, especially inductive processes 
neutralizing the total electric current, in the present study 
we use a general and fully self-consistent (basic plasma 
physics) approach - a 3-D electromagnetic particle-in-cell 
(PIC) modelling. 

All the abovcmcntioned processes necessarily modify 
the electron distribution function in the flare X-ray source. 
Moreover, contrary to simple models, which generally pre- 
dict high anisotropy of electrons and X-rays, it was found 
that the observed hard X-ray directivities are low (e.g. Kane 
1983). Furthermore, Kontar & Brown (2006) found a low 
anisotropy of the electron distribution function in the X- 
ray source by separating the reflected X-ray emission from 
the direct one. They concluded that the conventional solar 
flare models with downward beaming are excluded. 

In the present paper we want to demonstrate the impor- 
tance of the abovementioned processes on the evolution of 
the beam-plasma system with the return current. Our aim 
is to show their effects on the anisotropy of the electron dis- 
tribution function in this system and thus on the directivity 
of the corresponding X-ray emission. Using the 3-D electro- 
magnetic PIC model, for the first time in the study of X-ray 
directivity, we compute the evolution of the beam-plasma 
system with the return current depending on the magnetic 
field in the beam propagation direction. Then, assuming 
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Table 1. Model parameters. 



Model 


mi/mo 




Vh/c 




A 


16 


1/8 


0.666 


0.0 


B 


16 


1/8 


0.666 


0.1 


C 


16 


1/8 


0.666 


0.5 


D 


16 


1/8 


0.666 


0.7 


E 


16 


1/8 


0.666 


1.0 


F 


16 


1/8 


0.666 


1.3 


G 


1 


1/8 


0.666 


0.0 


H 


1 


1/8 


0.666 


1.3 


I 


100 


1/8 


0.666 


0.0 


J 


100 


1/8 


0.666 


1.3 


K 


16 


1/8 


0.333 


0.0 


L 


16 


1/8 


0.333 


1.3 


M 


16 


1/40 


0.666 


0.0 


N 


16 


1/8 


0.234^ 


0.0 





16 


1/8 


0.234~ 


1.3 



" mean velocity of the power-law beam distribution 
' mean velocity of the power-law beam distribution 



that the resulting electron distribution functions generate 
X-ray bremsstrahlung, we calculate the directivity of the 
associated X-ray emission. (For a detailed analysis of the in- 
stabilities and waves produced in the studied beam-plasma 
system, see Karlicky et al. 2008, Karlicky 2009, Karlicky 
and Barta 2009.) 

The layout of the paper is as follows: In Section 2 we out- 
line our model. The results of computations of the electron 
distribution functions with the return current are shown in 
Section 3. In Section 4 we present the corresponding X-ray 
directivities. Finally, in Section 5 the results are discussed 
and conclusions given. 



Table 2. The real spatial and time scales as a function of 
the chosen plasma density ric. 



Tic 




t = 200/a;pc 


Ad 




(cm-^) 




(s) 


(cm) 


(s) 


10* 


5.64 X 10** 


3.55 X 10"' 


3.19 


3.61 


10^ 


1.78 X 10^ 


1.12 X lO"'^ 


1.01 


0.36 


10" 


5.64 X 10^ 


3.55 X 10"® 


0.32 


0.03 


10" 


1.78 X 10"' 


1.12 X 10~® 


0.10 


0.003 



power-law distributions are derived as mean distributions 
over the whole X-ray source for much longer timescales than 
those considered in the present study. In much smaller flare 
volumes and on much shorter timescales, the monoener- 
getic beam is a reasonable choice. Nevertheless, in Models 
N and O we added computations with the beam having a 
power-law distribution function. To show effects of instabil- 
ities distinctly we chose its power-law index (in the velocity 
space) as 1.5, and the low-velocity cutoff of 0.09 c. 

To keep the total current zero in these models in the 
initial states, we shifted the background plasma electrons 
in the velocity space (i.e. we initiated the return current) 
according to the relation Vd — — Wbfib/^'o, where Wb is the 
velocity of the electron beam, rib and ric are the beam and 
background plasma densities (for this type of initiation see 
Niemiec et al. 2008). The beam velocity was chosen to be 
Vh/c = 0.666 or 0.333 (in the z direction), see Table 1. The 
ratio of the beam and plasma densities was taken as nb /n^ 
= 1/8 (Models A-L and N-0), and n\,/n^ = 1/40 (Model 
M). 

Because computations in the PIC models are dimen- 
sionless, the results are valid for a broad range of plasma 



2. Model 

For our study we used a 3-D (3 spatial and 3 velocity com- 
ponents) relativistic electromagnetic PIC code (Buneman 
1993). The system sizes are = 45 A, Ly = 45 A, and 
= 600 A (where A is the grid size). 

For a basic set of models we initiated a spatially homo- 
geneous electron-proton plasma with the proton-electron 
mass ratio m-p/mc—lQ (Models A-F, and K-0 in Table 1). 
This is unrealistic and it was chosen to shorten the pro- 
ton skin depth and computations. Nevertheless, the ratio 
is still sufficient to well separate the dynamics of electrons 
and protons. For comparison we added models with the 
mass ratio mp/me=l and 100 (Models G-J in Table 1). 
The electron thermal velocity is vtc — 0.06 c (the corre- 
sponding temperature is Te = 21.4 MK), where c is the 
speed of light. In all models, 160 electrons and 160 protons 
per cube grid were used. The plasma frequency is Wpe = 
0.05 and the electron Debye length is Ad = 0.6 A. In the 
models with the proton-electron mass ratio TOp/mo=16, the 
electron and proton skin depths are Acc = 10 A and Ad — 
40 A, respectively. 

Then, we included one monoenergetic beam homoge- 
neous throughout the numerical box (see Models A-M). 
Note that due to the physical and numerical simplicity 
and the propagation effect in which faster electrons escape 
from the slower ones, in most cases we consider monoener- 
getic electron beams, although in the interpretation of solar 
flare hard X-rays, the power-law distributions are used. The 
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Fig. 1. The electron distribution functions in Model B at 
four different times: at the initial state (a), at Wpci = 40 
(b), at ujpct = 100 (c), and ujpct = 200 (d). Crosses cor- 
respond to f{vz), dotted and dashed fines display f{vx) 
and f{vy), respectively. Note tfiat f{vx) and f{vy) overlap. 
Tfie single cross in tfie part a) at w/c = 0.666 denotes the 
monoenergetic electron beam. 
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densities. The real time and spatial scales are given by spec- 
ifying the plasma density. Table 2 summarizes temporal and 
spatial scales (the interval of computations t = 200/ti'pc and 
the Debye length) for the plasma densities in the 10*-10^^ 
cm~^ range. The processes under study are very fast. The 
coUisional processes are much longer, see the collisional free 
time {1/vq) in Table 2. The numerical system size is small 
(45 A X 45 A X 600 A = 75 Ad X 75 Ad x 1000 Ad, i.e. for 
the plasma density e.g. Ue = 10^ cm~'^ it gives 76 cm x 76 
cm X 1010 cm). Since the periodic boundary conditions are 
used, in reality the studied problem is infinite in space. 

The beam density and the corresponding beam energy 
flux is given by the chosen plasma density rio, rib/n-e = 
(1/8 and 1/40), and the beam velocities (see Table 1). For 
example, for rif. = 10^ cm~'^, nb/ug ~ 1/8, and Vf, = 0.666 
c, the beam density rib — 1-25 x 10^ cm^^ and the beam 
energy flux E'flux — 4.55 x 10^^ ergs s~^ cm~'^. 

Because we want to study the influence of the magnetic 
field, in the models we consider several values of the ratio 
of the electron-cyclotron and electron-plasma frequencies 
(wce/wpo = 0.0, 0.1, 0.5, 0.7, 1.0, and 1.3 - see Table 1). 
Note that in the space close to the flare acceleration site 
in the low corona there is plasma of relatively low density. 
Thus, for the huge electron beam fluxes required for an 
explanation of the observed X-ray bursts, such high ratios 
of nb/fio are needed. In all models, the periodic boundary 
conditions were used. 

3. Results of 3-D PIC simulations 

As an illustration of the time evolution of the electron dis- 
tribution function in the beam-plasma system with the 
return current. Fig. 1 shows this evolution for Model 
B. As can be seen, due to the two-stream instability 
(Michailovskij 1975), a plateau of the distribution func- 
tion f{vz) (in the beam propagation direction) on the beam 



side is formed. Moreover, some small part of the electrons 
even increased their energy due to their interaction with 
generated Langmuir waves. Simultaneously, the distribu- 
tion functions f{vx) and f{vy), i.e. the distribution func- 
tions in the directions perpendicular to that of the beam 
propagation, are strongly heated. This is due to the Wcibcl 
instability (1959) (see also Nishikawa at al. 2006). 

To demonstrate how the magnetic field influences the 
resulting electron distribution function. Fig. 2 presents 
the distribution functions for six values of the ratio 
of the electron-cyclotron and electron-plasma frequencies 
(tjco/wpo = 0.0, 0.1,0.5, 0.7, 1.0, and 1.3 - Models A-F, Table 
1). It is evident that with the increase of the ratio Wco/wpc, 
the role of the Weibel instability is more and more reduced, 
the distribution functions in the direction perpendicular to 
the beam propagation f{vx) and f{vy) are less heated. On 
the other hand, the problem of the return current forma- 
tion becomes more and more one-dimensional and a more 
extended tail on the return current side is formed (com- 
pare Model A and F in Fig. 2, see also Karlicky et al. 2008; 
Karlicky 2009). In Fig. 3 the same results are expressed in 
terms of the electron distribution functions depending on 
the electron energies. Although this type of description is 
more common in flare research, the distribution functions 
in velocity space presented in Fig. 2 carry more information 
than those in Fig. 3 and thus they are more physically rel- 
evant in describing the studied processes. The ratio of the 
electron kinetic energies in the direction parallel and per- 
pendicular to that of beam propagation, which expresses 
the "anisotropy" of the system, is shown in Fig. 4. The 
ratio of energies is defined as: 

where n is the number of electrons in the whole numeri- 
cal box. As can be seen in Fig. 4, the coUisionless (wave- 
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Fig. 3. The electron distribution functions in electron energies (thick lines) at iVpet = 200 as a function of the magnetic 
field in Models A-F with Wce/wpe = 0.0, 0.1, 0.5, 0.7, 1.0, and 1.3, respectively. For comparison in each panel the initial 
electron plasma distribution is added (thinner lines). 
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Fig. 4. Time evolution of the ratio of the electron kinetic parallel and perpendicular energies Epa,r/Epei as a function of 
the magnetic field in Models A-F with Uce/cOpe = 0.0, 0.1, 0.5, 0.7, 1.0, and 1.3, respectively. 



particle) processes very rapidly decrease the "anisotropy" 
on time scales shorter than Wpot ~ 50. This process is faster 
and more efficient for lower magnetic fields. While the end- 
ing ratio is i?par/£^perp ~ 9 for Model F (wcc/t^po = 1-3), in 
Model A (wce/wpe = 0.0) this ratio is only -Epar/-Eperp ~ 2. 



In Fig. 5 a comparison of models with three different 
mass ratios (mp/me — 1, 16, 100) and two values of the 
ratio oJcc/f-^pc (0.0 and 1.3) is made. While in the cases 
with TOp/rrie = 1 (the electron-positron plasma) the strong 
heating of the distribution functions /(f^) and f{vy) can be 
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Fig. 5. The electron distribution functions at uj^J. = 200 as 
a function of the mass ratio: mi/m-e = 1 - two upper plots, 
rrii/mc = 16 - two middle plots, and m{/me = 100 ~ two 
bottom plots for two values of uOcc/^-pc = 0.0 (left column) 
and 1.3 (right column). Notation is the same as in Fig.[TJ 



Fig. 6. The electron distribution functions in Model N with 
the power-law beam at four different times: at the initial 
state (a), at oj-pct — 40 (b), at Wpci = 100 (c), and Lo-pct — 
200 (d). Notation is the same as in Fig.[TJ 



seen even for the strong magnetic field (cjco/wpo = 1.3), for 
the proton-electron plasma the resulting f{vx) and f{vy) 
for nip/me = 16 and 100 do not differ significantly. Note 
that in the model with mp/rric — 100 the proton skin depth 
is greater than the system sizes and Ly. 

We also compared the evolution of the electron distri- 
bution functions in Models A and F with Models K and L, 
i.e. the models with a lower initial beam velocity {vh/c = 
0.333). We found that only the extent of the return-current 
tail in Model L is shorter than that in Model F. It is a nat- 
ural consequence of the greater beam velocity in Model F 
than in Model L. Furthermore, it was found that Model M 
gave qualitatively the same results as Model A. 

In Figs. 6 and 7 the electron distribution functions in 
Models N and O, i.e. in the models with the power-law 
beam and with two different ratio of electron-cyclotron 
and electron-plasma frequencies (wce/wpe = 0.0 and 1.3) 
are shown. Because these models are not subject to the 
bump-on-tail instability there are no significant changes in 
the distribution J{vz) on the beam distribution side. On 
the other hand, the Wcibel instability plays its role, espe- 
cially in the case without the magnetic field (Model N). 
Once again, in Model N the plasma is heated in the direc- 
tion perpendicular to that of beam propagation, whereas 
in Model O, the return current is formed by the extended 
distribution tail. 
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Fig. 7. The electron distribution functions in Model O with 
the power-law beam at four different times: at the initial 
state (a), at uo^ct — 40 (b), at Wpci = 100 (c), and ujpct = 
200 (d). Notation is the same as in Fig.[TJ 



4. Directivity of X-ray emission 

Knowing the electron distribution function /(v), an in- 
stantaneous X-ray bremsstrahlung, i.e. the so-called thin- 
target emission (e.g. Brown et al., 2003) can be calculated. 
To account for the anisotropy of f{v), we considered the 
angle-dependent electron-ion bremsstrahlung cross-section 
(5(e, E, 8) differential in the electron energy E and the solid 



angle of the incoming electron, where e is the photon en- 
ergy and 8 is the angle between the electron pre-coUision 
velocity and direction of the photon emission (Gluckstern 
& Hull, 1953). We used the expression for Q(e,E,9) given 
in Appendix of Massone et al. (2004), which includes the 
Elwert (1939) Coulomb correction. The cross-section was 
evaluated using hsi_reg_ge_angle_cross . pro available in 
the Solar Software. 




Fig. 8. The X-ray directivity in several energies for f{v) corresponding to Models A and F at ujpct = 200 and the X-ray 
directivity in the initial state for Models A-F (the case of simple beaming). The horizontal solid line represents the 
isotropic case, the dashed vertical line denotes the viewing angle for a limb source. 
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Fig. 9. The electron directivity in several energies for Models A and F at uipct = 200. The horizontal solid line repre- 
sents the isotropic case, the dashed vertical line denotes the viewing angle for a limb source. The corresponding X-ray 
directivities are shown in Fig. [51 



Figure [5] shows the X-ray directivity, i.e. the ratio of 
the angle-dependent J(e, 9) to integral photon spectrum 
/(e) = 1/Att f^I{e,9,(j)) dfi, where 9 and is the polar 
and azimuthal angle, respectively, is the solid angle. The 
z-axis of the coordinate system is chosen to be along the 
beam propagation direction. Note that due to axial symme- 
try of the problem around the z-axis, the photon spectrum 
I{e,9,(j)) is also independent of </>, so I{e,9) = I{e,9,(f). 
Assuming that the beam propagates along the local normal 
line towards the photosphere, Fig. [5] displays a variation of 
the X-ray directivity observed from different viewing an- 
gles: the cases with cos 9 — 1 and cos 6' = — 1 correspond to 
the forward (the direction to the photosphere) and back- 
ward (the direction to the Earth's observer when the X-ray 
source is at the disc centre) emissions, while the case with 
cos = denotes the emission in the perpendicular direc- 
tion (the X-ray source placed on the solar limb). 

The behaviour of the X-ray directivity is closely re- 
lated to the corresponding electron distribution. Comparing 
Model A and F at the time uj^^t = 200 with Models A-F 



in the initial state (i.e. the case with a simple beaming) 
in Fig. [51 it can be seen that values of the directivity, es- 
pecially in the backward direction, become closer to the 
value 1 (the isotropic case). Therefore, the global directiv- 
ity decreased during the evolution of the electron distribu- 
tion. Furthermore, we can see that the directivity values 
for cos 6* = in Model A are closer to the isotropic case 
than those in Model F. This is due to the strong heating 
of the plasma in the direction perpendicular to the beam 
propagation and it is caused by the Weibel instability in 
Model A (the case with zero magnetic field). 

We also defined the electron directivity f{E,9)/f{E), 
similarly to the X-ray one. Models A and F at time uipct 
— 200 are presented in Fig. [9l and show in another way 
the electron distribution characteristics discussed above in 
Section 3, Fig. [2] Comparing these electron directivities, we 
can see that they differ more distinctly than the correspond- 
ing X-ray directivities (Fig. [8]). Such a difference is caused 
by the strong smoothing effect of the bremsstrahlung cross- 
section. 
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We also calculated the X-ray directivities for Models K- 
L and N-0. They show the same changes as follows from 
the comparison of plots in Fig. 8, but these changes are less 
pronounced due to smaller changes of the f{v) anisotropy 
in these models. 

5. Discussion and conclusions 

Varying the ratio of electron-cyclotron and electron-plasma 
frequencies Wcc/'^pc, it was found that the magnetic field in- 
fluences the evolution of the electron distribution function 
in electron beam - plasma system with a return current. 
While for small magnetic fields (wco/wpo < 0.1) the elec- 
tron distribution function becomes broad in the direction 
perpendicular to the beam propagation due to the Weibel 
instability and the return current is formed by the elec- 
trons in a broad and shifted bulk of the distribution, for 
stronger magnetic fields (wco/^^pc > 1) the distribution is 
more extended in the beam-propagation direction and the 
return current is formed by the electrons in an extended 
distribution tail. Assuming the magnetic field and electron 
density a.s B = 100 G and no = 10^^ cm~'^ relevant to so- 
lar flares, the ratio of the electron-cyclotron and electron- 
plasma frequencies is Wce/wpe = 0.1. In such conditions the 
Weibel instability plays a role, but it is reduced for a higher 
magnetic field. The evolution is influenced also by the two- 
stream instability. Besides the formation of the plateau of 
the electron distribution on the electron beam side, the si- 
multaneously generated Langmuir waves even accelerate a 
small part of the electrons. 

The coUisionless processes cause a very fast decrease of 
the ratio of the electron kinetic parallel and perpendicu- 
lar (with respect to the beam propagation direction) en- 
ergies and lead to a decrease of the "anisotropy" of the 
system. Thus, the distribution function rapidly deviates 
from that with simple beaming. This can be also expressed 
by a decrease of the directivity of the associated X-ray 
bremsstrahlung emission. This fact agrees with the state- 
ment of Kontar & Brown (2006) that conventional solar 
flare models with a simple downward beaming should be 
excluded. 

An additional aspect of the present study is that the 
inclusion and physical necessity of the return current in 
the beam ~ plasma system resolves the problem of num- 
ber of electrons needed for an acceleration of the dense 
electron beam in the corona where the density is relatively 
low. The return current simply carries the same amount of 
electrons as in the electron beam back to the acceleration 
site. However, the return current does not have the same 
distribution function as the initially injected beam. 

Variations of the X-ray directivity obtained in our mod- 
els are of a level comparable to those in the electron beam 
propagation models by Langer & Petrosian (1977, Fig. 1) 
and Leach & Petrosian (1983, Fig. 4). However, there is an 
important difference betwee n o ur model and the models by 



Langer & Petrosian (1977) and Leach fc Petrosian (1983) 
We treat only coUisionless processes which were neglectec 
in the previous studies. Due to the very short time scales 
in our computations, no effects of longer beam propagation 
or collision scattering are included in the electron beam 
evolution. 

Therefore, the similar level of X-ray directivies suggests 
that a comparable level of isotropisation of the electron dis- 
tribution function caused by the coUisional processes can be 



produced by the studied wave-particle processes on much 
shorter time scales. Moreover, it means that these fast pro- 
cesses should not be neglected in X-ray directivity studies. 

Our study is not aimed at a direct comparison with 
observations, mainly due to the large difference be- 
tween simulated and observationally available time scales. 
Nevertheless, the paper by Kontar & Brown (2006) allows 
us to compare our simulations with their derived ratio of 
downward-to- upward electron distributions, Fd{E) / Fu{E) . 
The comparison reveals an agreement between inferred 
Fa{E)/F^{E) and Model F within the confidence interval 
up to ^ 50 keV. At higher energies, our models predict a 
directivity higher than that obtained from observations. 

The results presented here could be appropriate for low- 
density parts of flare loops where the coUisionless processes 
are dominant. Furthermore, one may consider them as in- 
put into simulations (on much longer time scales) which 
treat a propagation of the beam in the environment where 
Coulomb collisions play a significant role, such as the tran- 
sition region and the chromosphere. Since all these pro- 
cesses (coUisionless on long time scales, coUisional and even 
ionization processes in the background plasma) lead to fur- 
ther isotropisation of the particle distribution, we speculate 
that the resulting electron distribution and X-ray directiv- 
ity would be much closer to the isotropic case, as was re- 
cently found from X-ray observations (Kontar & Brown, 
2006). 
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